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E'JDLEY K.'rox LIEH'.,- 
^AVAL rr, ■ - , ■ 

~ -Y. C,.L-F.' ’ I 

ABSTRACT 

Three methodologies from which to approach the problem 

of preliminary helicopter design are explored in this paper 

The first is a sensitivity analysis of the basic helicopter 

performance equations. The purpose here is to ascertain 

where reasonable simplifications can be made that do not 

seriously degrade the accuracy of the results. The second 

is a graphical parametric design method, known as Carpet 

Plots. In this method a graphical solution is developed 

to meet the design criteria of the helicopter. In the 

third, an overview of Boeing Vertol's Helicopter Sizing and 

Performance Computer Program is given. The computer 
• • 
routines which enable a person to access HESCOMP on the 

Naval Postgraduate School main frame IBM system are also 

provided . 
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I. INTRODUCTION 



A. GENERAL 

The helicopter design process, the subject of numerous 
articles and studies is an evolving discipline that borders 
on being an art. A successful design must balance the 
user's needs and desires against practical capabilities. 

With the introduction of composite materials and new 
technologies, principally in rotor and engine performance, 
significant advances have been made in helicopter capabil- 
ities. In some instances, the performances of hybrid 
helicopter designs rivals that of a similarly sized conven- 
tional aircraft. For example, the YVX, a joint Boeing-Bell 
venture, will have the hover and low speed capabilities 
of a helicopter while being able to cruise at 300 knots. 

Viable commercial and military helicopter designs are 
only thirty years old. The first major use of helicopters 
occurred during the Korean conflict. To put this in 
perspective, the first large scale use of conventional type 
aircraft was in World War I. 

Helicopter design can proceed on a number of different 
levels, ranging from comprehensive computer design programs 
to preliminary analysis using simplifications of the basic 
performance equations. Each has its merit and place. 
Computer-aided design provides a great deal of data. 
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Generally, these programs integrate aircraft configuration 
sizing, performance and weight calculations in an iterative 
process. An example of a computer design program for heli- 
copters is the Helicopter Sizing and Performance Computer 
Program [HESCOMP] , orginally developed by Boeing-Vertol 
for NASA. This program is currently used as a wide number 
of institutions conducting studies in helicopter design. 

On the opposite end of the spectrum would be sensitivity 
design studies using the performance equations. Surprisingly 
accurate simplications of these equations can be made. This 
provides the designer with an excellent method for doing 
first cut preliminary helicopter sizing at a low cost. 

B. OBJECTIVE 

This report is an investigation of several of the 
methods employed in the preliminary design of a helicopter. 
Conceptually, the report can be divided into three parts. 

In the first section, a sensitivity analysis of the basic 
performance equations is performed. The purpose here is to 
ascertain where reasonable simplications can be made 
that do not seriously degrade the accuracy of the result. 

In the second section a graphical method of doing 
parametric design studies, known as Carpet Plots, is 
developed. This method allows the user to formulate a 
graphical solution matrix to meet the design criteria 
specified for the helicopter. Carpet Plots are 
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particularly instructive since they give visual insight 
into the interplay of the various design parameters. 

In the last section, an overview of HESCOMP is given. 
Programs are developed which enable a person to access 
HESCOMP on the Naval Postgraduate School Main Frame IBM 
system . 
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II . SENSITIVITY AxMALYSES OF BASIC HELICOPTER EQUATIONS 



A. DESCRIPTION OF PROBLEM 

In preliminary helicopter design, there are a number of 
instances where a quick first cut analysis would be 
extremely helpful. This is especially true in determining 
the preliminary size of the helicopter required to meet 
the specifications. 

Historically, there are a number of variables in the 
performance equations of helicopters which may be treated 
as constants. This may allow for significant simplifica- 
tions and aid in the preliminary design process. 

In this section, a sensitivity analysis of the per- 
formance equations is done. In a sensitivity analysis, 
each parameter [or variable] is varied in order to deter- 
mine its effect on the equation. Variables which are 
shown to have little effect may be treated as constants and 
the equation simplified accordingly. 

B. SOLIDITY 

Solidity, a , is the fraction of the disk area that is 
composed of blades. It is a function of b , the number 
of blades, of a constant cord, c , at a radius, R: 
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C. DISK LOADING: 



Disk loading is defined as the ratio of the weight 
the total area of the rotor disk. 

nr _ WEIGHT 
AREA 



to 



( 2 . 2 ) 



W 

A 




[Ib/ft^] 



D. POWER LOADING 

Power loading is the ratio of weight to input power. 



PL = [Ib/hp] (2.3) 

in 

In a hover, thrust equals weight; this allows us to 
rewrite the power loading for the hover condition as 



PL 



T _ ROTOR THRUST , 

P. ROTOR HORSEPOWER L-LD/npj 
in 



(2.4) 



E. COEFFICIENT OF THRUST AND POWER 

The coefficient of thrust, C,p , is a non-dimensional 
coefficient which facilitates computations and comparisons: 



T ^ T 

ApV.j.^ ttR^p(^^R)^ 



(2.5) 



Similarly, a coefficient of power, 
established as : 




has been 
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p 



p 



( 2 . 6 ) 



P ApV^-^ 






No significant simplifications can be made to either of 
these coefficients. However, it should be observed that 
the coefficient of thrust is inversely proportional to the 
square of the rotor tip velocity, while the coefficient of 
power is inversely proportional to the cube. 

Assuming all other factors being equal, increasing the 
rotor tip velocity from 600 fps to 700 fps [an increase of 
16.7 percent] will have the following result on these 
coefficients. 




T 

Ap(1.167)" 

T 

Ap(1.361) 



(2.5) 



The coefficient of thrust is reduced by 26.9 percent. 
Similarly, for the coefficient of power: 
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p 



c 



p 



3 



ApV^ 



P 



( 2 . 6 ) 



Ap(1.167)^ 



P 



Ap(l. 589) 



The coefficient of power is reduced by 37.1 percent. 
F. HOVER POWER 

The total power in a hover is made up of two terms, 
profile power, P^ , and induced power, P^ . 

Utilizing black element theory the prof ile , power 
required to hover can be expressed as: 



The induced power predicted by momentum theory is: 






(2.7) 




( 2 . 8 ) 




/' 



2it p R 



2 



The total power required to hover is: 




(2.9) 



16 



p, 



T 




( 2 . 10 ) 



Donald M. Layton in Helicopter Performance , [Ref. 1], 
found that for the optimum hover power, the induced power 
is equal to twice the profile power. The analysis was 
performed in the following manner. 

By assuming constant weight, density, solidity, and an 
average profile drag coefficient, as well as a fixed 
rotational velocity, equation (2.10) reduces to 



where and C 2 are constants. 

As equation (2.12) shows, profile power increases as 
the square of the blade radius while the induced power 
decreases with increasing blade radius. 

The optimum hover power with respect to rotor radius 
can be determined by taking the differential and setting it 
equal to zero. 




( 2 . 11 ) 





(2 .12A) 



or 




2 



(2.12B) 



which implies 



P^ = 2 P 



o 



( 2 .120 
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G. HELICOPTER SIZING 



A simplified relationship between the total power 
required, gross weight and rotor radius can be developed in 
the following manner. 

The total power required to hover equation for the main 
rotor was developed in the preceding section and is 
repeated here for clarity. 



Pt = P- + P 
T 1 o 



(2.9) 



Pt = 



^ 3/2 

/2tt p 



1,1_P .. 3„2 

R * 8 °r ‘^do f’' Aip ^ 



( 2 . 10 ) 



In a hover, thrust equals weight. Solving equation 
(2.11) for weight one obtains: 

= [P^ - i a pTT R^] (2.13) 



This equation may be further simplified if it is 
assumed that the density, average profile drag coefficient 
and tip velocity are constants; these are reasonable assump- 
tions. Historically, the average profile drag coefficient 
of a helicopter has been approximately 0.01. The operating 
environment of today's helicopters, especially military, 
is below 5,000 feet agl. This allows for the use of the 
standard sea level value for density with little error. 
Primarily, due to tip mach effects, the upper limit on the 
rotor tip velocity is in the range of 700 fps . 
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The resulting equation v\^ith these assumptions incor- 
porated into a constant, K , is: 

W = [47.527 P.J. R - bc]^/^ (2.13) 



Equation 
o£ magnitude 



(2.13) can be further reduced when the order 
of the two terms is considered. 



47.527 R >> be 



Thus , 



W [47.527 R]^/^ (2.14) 

To determine how accurate this simplification is, the 
equation is used to approximate the total weight of a 
number of helicopters for which the parameters are available. 
As Table 2.1 indicates, the weight approximation formula 
yields values within six percent of the actual total weight 
of these helicopters. 

H. FIGURE OF MERIT 

A figure of merit, FM , has been defined for the 
helicopter as the ratio of the ideal rotor induced power to 
the actual power required to hover, with non-uniform induced 
velocity, tip losses and profile drag power. 
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TABLE 2.1 



HELICOPTER WEIGHT COMPARISON 



HELICOPTER 


TOTAL GROSS 
WEIGHT 
(1000 lbs) 


CALCULATED 
GROSS WEIGHT 
(1000 lbs) 


PERCENT OF 
ACTUAL 

GROSS WEIGHT 


AH-64 


14.66 


14.69 


1011 


UH-IN 


14.20 


13.74 


97% 


H-3H 


21.00 


20 . 63 


98% 


S76 


10.00 


9.90 


99% 


UH-6DA 


20.25 


19.33 


95% 


H-543 


42.00 


42.00 


100% 


H-53D 


42.00 


41.00 


98% 


H-53E 


73.50 


69.00 


84% 
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In a hover, the figure of merit may be written as: 



FM 



1 

s/T~ 



. PL . 



DL 

/~ 



CT^ • ^ 
/2C 

P 



(2.15) 



The figure of merit is customarilty plotted against the 
quantity CT/a . According to Zalesch [Ref. 2], CT/a , is 
proportional to the average blade angle of attack and can 
be used as a measure of rotor efficiency. The curve in 
Figure 2.1 is based on data from Reference 2 for a typical 
tail rotor helicopter. 

Main Rotor Hover Performance 




Figure 2.1. FM Versus Blade Loading CT/a 
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Previous studies have shown that a figure of merit 
between 0.70 and 0.80 is considered average. [Ref. 5] 

If the induced power is between 70 and 80 percent of the 
total power, the figure of merit will be approximately 0.75. 

With the figure of merit limited to values between 
0.70 and 0.80, the following simplification can be made, 
assuming the hover condition of thrust equaling weight and 
standard sea level conditions; 



FM 



W 



3/2 



67.214 P.j,R 



(2.16) 



For Navy helicopter design, the rotor radius has been 
limited by flight deck spotting constraints to less than 
30 feet; the exception to this is the H-3, R = 31 feet and 
the H-53, R = 36 to 38 feet [depending on the model]. 
However, these two helicopters work almost exclusively from 
large air dedicated ships such as the LPH, LHA and CV. 

If the small deck operating assumption is made, 
equation (2.16) can be further simplified to [assuming 
R = 28 feet] : 

W ^ ^ 

^ " 1881.98 FM (2-17) 

An FM of 0.80 will yield a P to W relationship of: 

1 ^ 3 / 2 

" 1505.58 (2.18) 
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v>:hile an M o£ 0.70 yields a relationship 



P 



„ 3/2 

1517.39 



(2.19) 



I£ equation (2.17) is solved utilizing the approximate 
weight relationship developed earlier o£ 

= 47.527 P^R (2. 14) 

a value £or the £igure o£ merit o£ 0.707 is obtained. This 
is within the historical range o£ values. 

I. TAIL ROTOR SIZING 

A historical analysis o£ typical helicopters [Re£. 3), 
shows the £ollowing empirical relationship £or the tail 
rotor radius 



^ 1.3 [£t] (2.20) 

when comparing the results o£ this equation with actual 
tail rotor radius data, it was £ound that i£ a multipli- 
cation £actor o£ 1.2 is used vice 1.5 a better approximation 
is obtained. The results are tabulated in Table 2.2. 



J. FORWARD FLIGHT POWER CONS I DEFLATIONS 

The total power in Forward Flight consists oF induced, 
proFile and parasite power. IF the helicopter is a single 
rotor vehicle, the tail rotor power should be taken into 
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TABLE 2.2 



TAIL ROTOR SIZING 



HELICOPTER ACTUAL TAIL ROTOR APPROXIMATION [FT] 





RADIUS [FT] 


[2.20] 


[2.21] 


AH-64 


4.6 


4.98 


4.59 


UH-IN 


4.3 


4.90 


4.52 


SH-3H 


5.3 


5.95 


5.5 


S-76 


4.0 


4.11 


3.79 


UH-60A 


5.5 


5.85 


5.4 


CH-55D 


8.0 


8.42 


7.78 


CH-53E 


10.0 


11.15 


10.29 
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account, as well as all mechanical losses [transmission, 
etc.] for accurate calculations. Hoi^rever, a reasonable 
approximation can be obtained by considering only the main 
rotor and increasing this power figure by several percent 
to account for these losses. 

Figure 2.2 is a plot of the induced, profile, parasite 
and total power curves for typical tail rotor helicopter. 

MAIN ROTOR POWER 



SSL CONDITIONS 
OUT OF GROUND EFFECT 




Figure 2.2. Pov\rer Required Versus Forward Velocity 

The induced power drops off rapidly with increasing 
forward- velocity , whereas the parasite power increases 
rapidly . 

Parasite power is the power required to overcome the 
drag forces created by the aircraft's geometry. These drag 
forces are due to pressure drag and skin friction. 
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Parasite drag is extremely sensitive to the helicopter's 
loading. It is generally a minimum for forward flight and 
increases for side;\rays flight. Helicopters are generally 
streamlined for forward flight and the flat plate area is 
a minimum in this direction. The equation for the parasite 
power is: 




The parasite power is a function of the cube of the 
forward velocity. As such, with the advent of high speed 
helicopters a great deal of consideration has been placed 
on streamlining the geometric shape in order to reduce this 
power requirement. 

Blade element theory is commonly used to develop the 
profile power equation for forward flight. An excellent 
development of this equation is given in Reference 1. 

The profile power equation in forward flight is: 

Pof = ^ a p A [1 + 4.3 p^] (2.22) 

Equation (2.23) is a function primarily of the main 
rotor geometry. The variable with the most significance 
is the rotor tip velocity; increasing the tip velocity from 
600 to 700 fps results in a 58.8 percent increase in 
profile potver [assuming other factors are constant]. 
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The induced power is a function of the induced velocity. 
In a hover, the total flow through the rotor system is 
induced. As the forvv'ard velocity increases, the mass flow 
rate through the rotor disc increases due to the forward 
translation of the helicopter. This reduces the induced 
velocity . 

The equation for the induced power requirements at all 
forward velocities is: 



P = T . V-^ 

It 



(2.23) 



where 



^it 






2 2 

/2V 



2 * 



1/2 



+ 1( .V (2.23a) 



At high forward velocities, the induced power required 
can be approximated as: 



W 

p. = wv.. = 

1 it 2pAV£ 



(2.24) 



The total power for forward flight is the sum of the 
induced, profile and parasite powers. 



Pt = Pi ^ Po ^ Pp 



(2.2 5) 
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Pt = T,v.^ ^ 1 ^ " 4.3 ,^3 



+ 2 P ff 



(2.25a) 



At high forv\/ard velocities, equation (2.23) can be 
substituted into equation (2.25), resulting in: 



2 ..V 

" * i * VT-qi * 4.3 



T 2pAV. 8 






* ^ P ff V/ 



If one makes the following assumptions 



W = const = const 



(2.26) 



p = const 
VT = const 



a = const 



Equation (2.26) reduces to 



2 2 

P„ = ^ R + 

T r2 p 



(2.27) 



The derivative of equation (2.27) with respect to radius is: 



dP.^ 

dF 



2K. 



^ . 2 K2 R 



(2.28) 
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Setting this equal to zero, one obtains: 



2X 



U W R = 0 



(2 . 28a) 



R 



R ^ 
2 



2K 



1 



+ 2 R] = 0 



(2.28b) 



K 




(2.28c) 



R 




(2 .28d) 



This defines point of minimum total power required for 
VMAX range. This corroborates with the res.ults obtained by 
Waldo Carmona [Ref. 4]. 

If the total power required is differentiated with 
respect to forward velocity and is set equal to zero, it can 
be seen that 



P. = 3 P 
1 



o 



(2.29) 



or 



w2 _ 3pfV^ 



2 



2pAV^ 2 



(2.30) 
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Solving this equation for velocity results in: 



V 



£ 




A 





ft/sec 



(2.31) 



According to Carmona [Ref. 4] , this corresponds to the 
best endurance velocity. 

K. DENSITY EFFECTS ON TOTAL POWER 

The effect of density on the total power required in 
for;\fard flight is as follows: 

The general operating altitudes of a helicopter are 
belo\\T 10,000 feet. The corresponding ICAO STANDARD 
ATMOSPHERE range for density is 

p = 0.0023769 [lb sec^/ft"^] SSL 
p = 0.0017553 [lb sec^/£t^] at 10,000 feet 

p/pSSL varies from 1 to .7385. 

The effect on the components of P.^. are as follows: 
Induced Power: 

1/p/pSSL => 1 to — 

This translates to a 35 percent increase in the induced 
power . 
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Parasite and Profile Power: 

Both parasite and profile powers are directly propor- 
tional to the density ratio. Therefore, as you go up in 
altitude both P 

o 



and P are reduced. 
P 
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Ill . CARPET PLOT DESIGN STUDY 



A. DESCRIPTION OF PROBLEM 

Preliminary helicopter design involves one with a wide 
range of choices. For any given payload and performance 
specifications, there a number of helicopter designs that 
satisfy the requirements. The problem in the preliminary 
design process is narrowing these possibilities and 
selecting the design which will provide the best helicopter 
for the mission. 

Obviously, the operating environmental constraints help 

to define the basic configuration. These constraints are 

usually specified in the Request for Proposal [RFP] , in the 
^ • 
case of a military helicopter. For example, typical 

constraints placed on the design of a Navy helicopter are 

the size of the ship deck and hangar from which it will be 

operating, the requirement for a blade fold system, dual 

engine configuration and IFR capability. 

Even with these design constraints, there is still 
a great deal of leeway. In order to insure that the best 
helicopter design is selected, an appropriate number of 
solutions satisfying the specifications should be inves- 
tigated. Since each solution is generally characterized 
by a different combination of design parameters, the 
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selection, according to Greenfield [Ref. 5] , can best be 
made through a parametric study which allows for the 
optimization of many design parameters. 

One method of parametric analysis used is Carpet Plots. 
This method is based on the simultaneous graphical solution 
of the weight and hover performance equations. To this 
solution set is added to the environmental constraints to 
the helicopters size. This effectively brackets the area 
of acceptable design solutions. 

This method assumes that minimum gross weight is the 
criterion by which the best [or optimum] design parameters 
are selected. 

B. ASSUMPTIONS 

1. Airfoil used is a derivative of the NACA 0012 

with the following mean approximate values from Reference 5. 
a = slope of airfoil section lift curve, dC^/da , 
per rad. 
a = 5.73 

S = blade section drag coefficient 
6 q = .009 

62= -3 

2. a) The tail rotor radius is assumed to be .16 
times the main rotor radius [Ref. 5]. 
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b) The distance betvvreen the rotors, or tail rotor 
moment arm, is 1.19R [Ref. 5]. These ratios reflect 

the values of maximum rotor diameter and overall length 
specified as size limitations. 

3. B = .97. Historical approximation [Ref. 7]. 

C . METHODOLOGY 

In order to properly develop the weight and performance 
equations required for a carpet plot design study, the 
payload and performance specifications of the helicopter 
are needed. This data is used to tailor the equations for 
the design. 

The equations will be developed here for a four-place 
light helicopter. The equation development procedure is 
applicable to other size helicopters; the development for a 
medium helicopter, 20,000 lb weight class, is to be found 
in Appendix B. 

The following specification requirements which are 
similar to those in Reference 5 will apply to this design: 

1. The rotor diameter should be less than 35.2 feet. 

2. The overall length should be less than 41.4 feet. 

3. The gross weight of the helicopter should not 
exceed 2,450 lbs. 

4. The helicopter should be capable of hovering, out 
of ground effect at 6,000 feet with an ambient air tem- 
perature of 95*^F. 
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5. The useful load at hover shall consist of, as a 
minimum, 200 lbs for the pilot, 400 lbs of payload and 
sufficient fuel to give the helicopter up to three hours 
endurance at sea level conditions. 

6. Maximum speed of at least 110 knots using Normal 
Rated Power, at sea level. 

7. Total Power Required at 6,000 feet and 95°F shall 
be not more than 206. 



D. HOVER EQUATIONS 

1. The main rotor power required to hover out of 
ground effect is 

Total iMain Rotor Power [Hover] = Rotor Profile Power + Rotor 
Induced Power 



1 . 1 5w rw~ 

550B/2^^ ./p/p^ 



6WV.J, 


p/Po 


. ^ [ ‘'LRo 


2- 


4400 


^LRo 


“o 1577^1 


. 



(3.1) 



At an altitude of 6,000 feet and a temperature of 95° , 
p/Po = .749395 . Therefore, equation (1) can be simplified 
to : 

Pt6000/95°F ■ •035479K[DL]1/2 

+ • 91 . 9 . 71 [lOj'^ri + 1.80779 C, ^ ^)W 
Ct LRo 1 

L Ro 
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The tail rotor thrust required to counterbalance the 



main rotor torque is: 



550 P, 



TR 



550 P^R 

1 . 19V^ 



T 



(3.3) 



where has been defined as 1.19R . With R^^^ defined 

as .16R , the tail rotor disk loading can be written, 
using equation (3) as: 



DL 



TR 



550 P. 



TR A, 



TR 



1.19V„ . ,...2 

T IT ( . 16R) 



(3.4) 



550 P. 



DL 



1. 19 ( .0256)V^ W 



Greenfield [Ref. 5], in his development, assumes that 
the tail rotor tip speed is equal to the main rotor tip 
speed and that 6.^^^ = .02 and = .90 . With these 

assumptions the equation for the tail rotor power required 
to hover can be written as: 

3/2 



= 2055.7 



TR. 



Hover 



DL 


1/2 


Pp 

Hover 


W p/pQ 




1 

> 

i 



(3.5) 



012605 P, 



Hover 



'LRTR 
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The equation for the tail rotor mean blade lift coeffi- 
cient can be written as 



C 



LRTR 562 . 5[p/p^) 



(3.6) 



if it is assumed that the tail rotor is designed to counter- 
balance a sea level main rotor torque equivalent to 90 
percent of the installed power. 

Substituting equation (3.6) into equation (3.5) one 
obtains the following expression for hover tail rotor povver: 



It is assumed that the gear losses amount to 3 percent 
and that there is a 1 percent cooling power loss, the total 
brake horsepower required to hover becomes: 

P + P 

^ _ ^Tm ‘^^TTR ,, 

96 

Empirical studies have shown that the tail rotor power 
required to hover can be approximated by 

P ~ .8 [total horsepower to hover] 

^AC 

This allows one to write the main rotor power required to 
hover as : 



P 



T, 



2374.7 ^ 



1/2 f P„ 



TT — + 5.3134 

Vt 



(3.7) 



TR6000/950 




(5.9) 
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Following Greenfield's [Ref. 5] development further, 
if equations (3.2) and (3.7) are substituted in equation 
(3.8) , one obtains 

^ = .0367 57 W /DT" 

^H6000/95° 



+ (10)’^ [1 + 1.80779 W V„ (3.10) 

Lrro LRo T 



+ 2473.6 





3/2 

+ 5.5348 



Utilizing the approximation for tail rotor power, 
equation (3.9), equation (3.10) can be solved for W 
(gross weight) as a function of variables V.p (tip speed) 



DL (rotor disk loading) , C 



LRo 



(rotor mean lift coefficient) 



and P„ (total power to hover) 



K. 



W = 



411. 51 



DL 



;/4 






1/2 






1 + K. 



/DL 



K. 



(3.11) 



/DL 



where : 



K-, = P 



( 10 ) 



1 ^16000/90" 



= ♦0 002 59 _ 2 _9 1.80779 



'LRo 



(3.11a) 



(3.11b) 
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K_ = 



5534S0 

K. 



(3.11c) 



K 



3695.7 

K. 



(3. lid) 



K, 



.95803 

^LRo 



(1 + 1.80779 C 



LRo 



(3. lie) 



Equation (3.11) has been programmed in Appendix B 
and solved for tip speeds from 600 to 700 cps and 
of .3 to . 7 . 

Equation (3.11) is one of the two primary equations used 
to obtain the data required for a carpet plot design 
analysis. Generally, the variables V.^, , DL , and 

P.J. , that are required for solution have specific ranges 
of values, depending on the weight class of the helicopter 
being designed. The graphical results of equation (3.11) 
for tip speeds of 600 to 700 fps and mean lift coefficients 
between .3 and .7 are illustrated in Figure 3.1. 

Both the Fortran and Disspla programs, as well as a 
decision making flow chart are provided in Appendix C 
to aid in using this method for a design solution. 

E. WEIGHT EQUATIONS 

Weight equations need to be developed that realistically 
reflect the sizing class of the helicopter being designed. 
The evolution is greatly simplified if a specific engine 
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Weight Equation Plot: CLR=0.5 
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Figure 3.1« Weight Fciuo^tion Plot, 



installation [# and horsepower] is assumed, since the 
weight of a number of components depend only on the 
installed power; this would include such terms as the engine 
controls and accessories. Another category would be those 
components whose weights depend on either the gross weight 
on two or more of the following in combination: rotor tip 
speed (V.J,) , rotor diameter (R) , rotor solidity (a) . 

The equations developed here are taken from the Hiller 
Aircraft Corporation Performance Data Report. [Ref. 5] In 
this report they assumed a specific engine installation, 
the Allison T-63 with a military power rating at sea level 
of 250 horsepower. 

There is a possible problem of the validity of these 
weight relationships when applied to different helicopter 
design categories. However, assuming a specific engine 
determines a number of the component weights, and thus 
minimizes the inaccuracies. Using the weight estimation 
relationships developed in the Helicopter Design Manual 
[Ref. 2], the engine, control and accessory weight can be 
calculated and the weight formulas developed here applied 
to give a representative useful load and empty weight 
formula for preliminary design analysis. This is done in 
Appendix C, for a 20,000 pound class helicopter. 
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The following relations are used to reduce the component 
weight formulas for the specification helicopter; 

W/DL = A = ttR^ (3.12) 



W/PL = MHP = 250 (3.13) 

(Military rating for Allison T-63 at sea level.) (PL = 

Power Loading.) 

P = (3.14) 



Using these equations the component weight for the 
specified helicopter empty weight may be reduced to the 
following : 



Engine, Controls 



Engine Section 
Group 



Main Trans- 
mission 



and Accessories 



.053 [W/PL] 



1.07 



10.43 



V/ 



1.295 



(PL V^) 



. 863 



= 617.5 lbs. 



= 19.5 lbs. (3.15) 

= 1221 p'^°^ (3.16) 



Rotor Drive ,^1.05 

Shaft 5.56 ^ 

(PL V^) ' (DL) 



266 p'"^ (3.17) 



Tail Rotor 



32.22 



W 



1.14 



(PL Yj) 



1.7 



17449 



V. 



1.14 



(3.18) 
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The engine, controls and accessories category includes 
such items as lubrication and oil cooling system, engines, 
communications, engine controls, engine accessories, 
instruments starting system, furnishing, flight controls, 
electrical system and stabilization. These are considered 
fixed weight items determined from specification of the 
engine and weight class of the helicopter. 



Tail Rotor 
Gear Box 



3.7 



W 



. 75 



(PL V.j.‘ ^(DL) 



59.47 /P 



(3.19) 



Tail Rotor 
Drive .124 

Shaft 



W 



1 . 355 



(PL V.p) ■ ^^(DL) • 



2.886 P‘^^ /A 



(3.20) 



Body and 

Gear 

Landing 



916 99 

= 1.91 W' + .0294 



Rotor ,,1.185 .33 ,, ,or -- 

Blade j5.15 ~ j5.15 A a 

Teetering Vy(DL)' ''T 



(3.21) 



Rotor Blade 
Artic- 
ulated 



,, ^^1. 205^. 33 ^ ,.205_.33 

19.77 - 19.77 TT — A o 

V„(DL)-"''^ '^T 



(3.22) 



Rotor Hub 71 

Teetering .0088 ~ -0088 WA' 

DL 



(3.23) 



Rotor Hub 21 

Artie- .00975 r-, = 00975 WA 



ulated 



DL 



. 21 



(3.24) 
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Fuel System .416 per gallon capacity = .0615 Wp (3.25) 

where = fuel weight. 

The individual component weights may now be combined 
into a single expression for the helicopter empty weight. 

W = 617.5 + .0617W^ = 1221P'^^^ + 266P'^ + il?-. 

e F ^1.14 

^ (3.26) 

+ 58.47/P + 2.886P'^^/A + .191W^^^ + .0294W^^ 

+ appropriate rotor blade and hub weights. 

As stated earlier, the design specifications called for 
a useful load consisting of a pilot (200 lbs), payload 
(400 lbs) and the required fuel weight (W ) . The fuel 

r 

weight is calculated for the Allison T-63 in the following 
manner: endurance of three hours at 85 percent of normal 
rated powered for the T-63 is 180.2 HP and the specific 
fuel consumption at this power is .783 lbs fuel/BHP HR. 
Including an allowance for a three-minute warm-up at NRP 
and using a 5 percent correction factor on SFC, as 
specified in Reference 5, the fuel weight becomes: 

Wp = 3(180.2) ( .822) + (212)(777) (3.27) 

.4n allowance should also be made for oil plus trapped 
fuel. This is estimated at 20 lbs. 
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The total useful load is the sum of the useful load 



items . 



= 200 + 400 + 452.6 + 20 = 1072.6 lbs (3.28) 

A nevv- variable, > is defined as the sum of the 

emply weight plus useful load. It is the of equations 
(3-2 6) and ( 3- 2 8) • 



Wear = 1717.9 + 1221P'^^^ = 266P'^ + + 58.47/P 



(3.29) 



57— 916 99 

+ 2.886P /A + .1911V + .0294W 



+ appropriate rotor blade and hub weights. 

Equation (3.11) together with equation (3.29) form the 
basis of a carpet plot design study. These equations are 
solved simultaneously for ’■'^gAR * This solution is best 
illustrated graphically, as in Figure 3.2. The graph in 
Figure 3. 2 was generated for a specific value of C^r over 
a range of tip speeds' [600 to 700]. 



F. GRAPHICAL .ANALYSIS 

Graphs similar to Figure 3.1 are generated for several 
value of C^R > and are then cross plotted to form Figure 

3.2. 



The mean lift coefficient, C^r > values are selected 
based on what is considered the historical average range of 
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Helicopter Carpet Plots: CLR=0.5 
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Disk Loading 

Figure 3.2. Helicopter Carpet Plots: 



values. Figure 3.5 is basic plot for a carpet plot design 
study. Programs are provided in Appendix D which will 
generate the required data sets and plots of Figures 5.2 
and 3.3. 

The solution field depicted in Figure 3.3 is too large 
to be of great analytic value and as such must be reduced. 
Three parameters, maximum gross weight, rotor diameter 
(both specified in the Design Specification) and the 
aspect ratio can be used to narrow the field of solutions. 

1 . Rotor Diameter Boundary 

A net to exceed value for the rotor diameter is 
generally given in the design specifications. This 
limiting value is based on the operating environment of the 
helicopter, h'ith R max specified, there is a linear 
relationship between the disk loading and the gross weight. 



BL = l = JL 
A p2 

TT R 



The resulting bracketing of the' solution field by 
applying both the maximum gross weight and maximum rotor 
diameter limits to the carpet plot are shown in Figure 3.4 



2 . Respect Ratio Boundary 

It is evident that a further restriction is still 
necessary to completely define the region of acceptable 
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Helicopter Carpet Plots 




(Sffi) W^ISM ssojg 



48 



Disk Loading 

Figure 3.3. Helicopter Carpet Plots 
Family of Solutions 



Helicopter Carpet Plots 




(sn) mSlsM SS0J3 



49 



Disk Loading 

Figure 3.4. Helicopter Carpet Plots Rotor 
Diameter and Weight Limits 



design solutions, 
rotor aspect ratio 
Thus 

21 
or 

For the case of a two bladed main rotor equation (3.30) 
reduces to: 

DL > .000012 

The detemination of this boundary graphically is as 
follows : 

2 

The hover solution plot of Figure 3.2 is replotted 
relative to the coordinates disk loading and design mean 
blade lift coefficient. The limiting curves for 
DL = .000012 are then plotted. The intersection 

with the appropriate constant tip speed lines of the hover 
solution represent the aspect ratio boundary; Figure 3.5. 



^For a helicopter rotor, the aspect ratio is defined as 
the radius divided by the chord. 

2 

For clarity lines of constant gross weight are omitted. 



Studies have indicated that a main 
of 21,^ is a representative upper limit 



R<mr> b ^ ^o *^LR ^T 



— C<mr> TTo 



OTT DL 



DL > 



b p 

126tt 
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aspect ratio boundary plot 
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Figure 3,5. Aspect Ratio Boundary Plot 



These intersection points are then cross plotted onto 
Figure 3.4. Figure 3.6 represents a graphical plot of the 
solution set satisfying the performance and structural 
design criteria of a small observation helicopter as 
specified in this study. 
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Disk Loading 

Figure 3 . 6 . Helicopter Carpet Plots 
Final Solution 



IV. 'HESCOMP 



A. DESCRIPTION OF PROGR.VM . 

HESCOMP is a helicopter sizing and performance computer 
program developed by the Boeing Vertol Company. The program 
was originally formulated to provide for rapid configuration 
design studies . 

A number of programming options are available to the 
user of HESCOMP. V/hen the type and mission profile of the 
helicopter are known, HESCOMP may be used to size the 
aircraft. Alternately, it may be used for mission profile 
calculations when the sizing details [gross weight, payload, 
engine size, etc.] are specified. A combination of these 
two options is also available; the program may be used to 
first size a helicopter for a primary mission and then 
calculate the off-design performance for other missions. 
Finally, HESCOMP may be used solely for obtaining helicopter 
weight . 

Sensitivity studies involving both design and per- 
formance tradeoffs can easily be done with HESCOMP. Incre- 
mental multiplicative and additive factors can be imbedded 
in the input data. 

The various helicopter configurations that may be 
studied using HESCOMP are detailed in Table 4.1. 
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TABLE 4.1 



HELICOPTER CONFIGURATIONS 
WHICH MAY BE STUDIED USING HESCOMP 



HELICOPTER COIJFIGURATIOnS HIUCII HAT BE STUDIED USING HESCOMP. 


t Iona 1 Ll f t/Propuln Ion 














tern Component F Whicli 














be AJJed to 




Propeller 




Type of Auxiliary 


Helicopter "Pure" 




for 


Auxl llary 


Independdnt Engines 


Type "'"'^-^Conf. 




Auxiliary 


Independent 
















(Doth Single ( Tandem Rotor) ^ 


Wing 


Propulsion 


Engines 


T/Shaft 


T/Fan 


T/Jet 


Pure Helicopter 






. 








Winged Helicopter 


X 












Compound Helicopter 














(1) Coupled (prim, engines 


X 


X 










drive auxiliary propulsion 
system) 














(2) Auxiliary Independent 














propulsion system 
(u) T/Sliaft engine 


X 


X 


X 


X 






(b) T/Fan cngitic 


X 




X 




X 




(c) T/Jet engine 


X 




X 






X 


Auxiliary Propulsion Helicopter 














(1) Coupled (prim, engines 




X 










drive auxiliary propulsion 
system) 














(2) Auxiliary independent 














propulsion system 
(a) T/Shaft engine 




X 


X 


X 






(b) T/Fan engine 






X 




X 




(c) T/Jet enqlne 






X 






X 


'oaxlai Rotor Helicopter 














(1) Coupled (prim, engines 




X 










drive auxiliary |iro|>ulslon 
sys tern) 














(2) Auxiliary Independent 














propulsion system 
(a) T/Shaft engine 




X 


X 


X 






(b) T/Fan engltie 






X 




X 




(c) T/Jet engine 






X 1 






X 
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B. progilajM modifications and implementation 



The computer program received from Boeing V''ertol 
required some modification and reformating in order to run 
properly on the Naval Postgraduate School IBM system. 

These alterations did not, however, alter the program output 
or usability. 

HESCOMP, as received from Boeing Vertol, uas 17821 
lines long and set-up as a sequential data set to be 
assemble on a 'G compiler'. The Batch processing system at 
the Naval Postgraduate School accepts only programs set to 
run on 'H compiler' . Normally, the differences between 
these two compilers are minor and programs that run on one 
will run on the other. However, this was not the case with 
HESCOMP . 

In order to facilitate the program debugging process, 
HESCOMP was reformatted as a partitioned data set. What 
this effectively did \vas to break the program down into 
eight members of approximately 2000 lines. The program 
breakdo;v’n is illustrated in Table 4.2. 

Each of these were compiled individually and then error 
codes analyzed. The member data set was then modified as 
required to properly compile. 

Once all the members of the partitioned data set 
compiled properly, HESCOMP was again formated as a 
sequential data set and run utilizing input data for 
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which there was a knovv^n output. This insured that the 
modifications made to the original program had not altered 
the logic, ie., gave faulty results. 

The control language program to access HESCOMB on the 
Batch processing system and a sample input and out data 
set are shown in Appendix D. These are also available 
on the Aero disk for copying and use. 

TABLE 4 . 2 



PARTITIONED DATA SET 



iER NAME 


LINE NUMBER 


SIZE 


FIRST ROUT] 


SI 


1 - 1681 


1681 


AERO 


S2 


1682 - 4132 


2451 


CLIMB 


S3 


4133 - 6531 


2399 


XIBIV 


S4 


6532 - 8974 


2443 


POWAVL 


S5 


8975 - 10870 


1896 


PRINT 1 


S6 


10871 - 13042 


2172 


ROT POW 


S7 


13043 - 15383 


2341 


CRUS 3 


S8 


15384 - 17821 


2448 


TAXI 


PROGRAM FLOW 


The program 


is conceptually 


outlined in 


Figure 4.1, 



[Ref. 7]. The program flow is monitored by a general loop, 
which controls a series of peripheral programs. There are 




Figure 4.1. HESCOMP Program Flow 
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a total of 44 subroutines. Detailed program descriptions 
cam be found in Section 4 of the HESCOMP User's Manual. 

D. PROGRAM INPUT 

Program input can be loosely group into ten categories : 
general information, aircraft descriptive information, 
mission profile information, rotor tip speed schedule, 
incremental rotor performance, auxiliary propulsion input 
schedule, engine cycle information, rotor performance 
information, propeller performance information, and 
supplementary input information. 

The actual amount of input data requires varies greatly 
with the program options selected. An example of a data 
set formatted to run on the IBM system is shown in 
Appendix E. A more detailed explantion is available in 
Section 5 of the HESCOMP User's Manual. 

E . PROGRAM OUTPUT 

An example of the program output is included in 
Appendix E. The printout consists of general data, input 
data, sizing data [program output] and mission performance 
data [for the size helicopter] . Detailed descriptions of 
these and diagnostic error statements are described in 
Section 6 of Reference 6. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



Three approaches to analyzing a preliminary helicopter 
design were explored in the course of this paper. It was 
found that a number of the performance equations could be 
greatly simplified with little degradation in the final 
results. A sensitivity analysis brought further insight 
into the inter-play of the parameters and how changes in 
them tended to effect the helicopter performance equations. 

Carpet Plots provided the most interesting method of 
analysis. Development of a graphical solution matrix 
using this method provides a usual interpretation of what 
is occurring when key parameters are varied. 

Two cases ivere explored; a light observation helicopter 
in the 3,000 pound weight class and a heavier utility 
helicopter in the 20,000 pound weight class. The Carpet 
Plot method provided reasonable solutions in both cases. 

In doing the analysis for the utility helicopter, the 
initial weight estimation equation had to be adjusted 
upward by approximately 2,000 pounds for the equations to 
intersect properly. This is not considered a limitation 
to this method of analysis, however, it does point up an 
area for further investigation. It may be possible to 
develop more accurate weighing factors for this equation 
when dealing with higher gross weight helicopters. 
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HESCOMP provides a plethora of information to the user. 
However, the price is the amount of inputed data required 
for even a simplified analysis. At a preliminary design 
level of analysis, the other methods explored provide a 
quicker first-cut look at the potential design. 
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APPENDIX A: NOMENCLATURE 



TERM 


DEFINITION 


UNITS 


a 


Slope of Airfoil Section 
Lift Curve 


Radians 


A 


Rotor Disk Area 


ft^ 


AR 


Aspect Ratio 


Dimensionless 


Atr 


Tail Rotor Disk Area 


ft^ 


b 


Number of Rotor Baldes 


Dimensionless 


B 


Tip Loss Factor 


Dimensionless 


C 


Main Rotor Cord 


ft 


n 

o 


Profile Drag Coefficient 
at Zero Lift 


Dimensionless 


^LRo 


Design Mean Blade Lift 
Coefficient at Sea Level 


Dimensionless 


Crj, 


Coefficient of Thrust 


Dimensionless 


C 

P 


Coefficient of Power 


Dimensionless 


5 


Blade Section Drag 
Coefficient 


Dimensionless 


DL 


Disk Loading 


Ib/ft^ 


FM 


Figure of Merit 


Dimensionless 


HP 


Horsepower 




Ltr 


Tail Rotor Moment Arm 


ft 


P 


Air Density 


lb sec^/ft^ 


P 


Advance Ratio 


Dimensionless 


R 


Rotor Radius 


ft 
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TERM 


DEFINITION 


UN I TS 


?! 


Total Power 


HP 


^TM 


Main Rotor Total Power 


HP 


p 

^TTR 


Tail Rotor Total Power 


Hp 


P 

0 


Profile Power 


HP 


P. 

1 


Induced Power 


HP 


P 


Parasite Power 


HP 


P 


PL 


Power Loading 


LB/HP 


R 


Rotor Radius 


ft 


T 


Thrust 


HP 




Induced Velocity 


ft/sec 


'^F 


Forward Velocity 


ft/sec 


V 


Aircraft Forward Speed 


ft/sec 


Vp 


Rotor Tip Speed 


ft/ sec 


W 


Aircraft Gross Weight 


lbs 


W 

c 


Empty Weight 


lbs 


“f 


Fuel Weight 


lbs 


w 

u 


Useful Load 


lbs 


‘'^BAR 


Empty Weight Plust 
Useful Load 


lbs 


a 


Solidity 


Dimensionless 
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APPENDIX B: CARPET PLOT FORMULATION FOR 20,000 LB. 

CLASS HELICOPTER 



B1 SPECIFICATIONS: 

Maximum Gross Weight: 20,000 pounds 

Maximum Rotor Diameter: 30 feet 

B2 PRELIMINARY ENGINE SIZING: 

B2 . 1 Utilize equation (2.14) to determine engine 
horsepower category. 



W = [4 . 753P.j,R] 
20,000 = [47.53P.J, 30]^^'^ 



P.J. = 1983 HP 

B2.2 Use the engine selection parameters tables B.l 
to determine the number and type of power plant 
[table taken from Reference 3] . 

B2.2a Type and number selected: 2 type C. 

B2.2b Specifications: 

Dry Weight Per Engine: 423 pounds 

Shaft Horsepower at Standard Sea Level: 

Military 1561 HP 

Normal 1318 HP 

B3 WEIGHT EQUATION FORMULATION 

B3.1 To obtain the engine control and accessory 
weight use items 7, 9, 10, 11, 12 and 13 of 
the weight estimation relationships developed 
in Reference 3 for a utility helicopter: 

#7: 609 lbs; “9: 129 lbs; #10: 76 lbs; 

#11: 410 lbs; #12: 439 lbs; and #13: 302 lbs. 
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TABLE B.l 



ENGINE SELECTION PARAMETERS 



The following turboshat power plant data are 
presented for one engine. 



Engines : 


A 


B 


C 


D* 


E 


F 


Dry Weight (lbs) 


158 


288 


423 


709 


580 


750 


SHP (ssl) Military 


420 


708 


1561 


1800 


2500 


3400 


Normal 


370 


659 


1318 


1530 


2200 


3000 


Cruise 


278 


494 


1989 


1148 


1650 


2250 


SFC (ssl) Military 


.650 


. 573 


.460 


. 595 


.615 


. 543 


Normal 


.651 


. 573 


.470 


.606 


.622 


. 562 


Cruise 


. 709 


. 599 


.510 


.661 


.678 


.610 


Initial Costs 


$93K 


SIOOK 


$ 580K 


$360K 


$640K 


$ 700K 


Operating Cost 
per hour/engine 


$8 


$16 


$20 


$35 


$40 


$60 


Preventative Maint 
per hour/engine 


$25 


$50 


$100 


$125 


$160 


$220 


MTBMA (hrs) 


3.5 


3.0 


2.0 


3.0 


4.0 


3.5 


MDT (hrs) 


0 . 7 


0.6 


0.5 


1.3 


2.0 


2.6 


MTBF (hrs) 


185 


210 


205 


285 


280 


320 


MTBR (hrs) 


600 


750 


800 


800 


1000 


750 
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B3.2 Simplifications 



B5 . 5 



B3.4 



B3. 5 



B3 . 6 



7 W "3 

k ■ ^ .■ iIe- "HP ■ 51.00 ; P - ^ ^ 

Engine Group 
.053(5100)^-°^ = 272 lbs 
Main Transmission 



10.43 



W 



1.295 



(^pm V^) 



. 8 b 3 , . .432 



= 10.43 



^^.863 ^+.432 
(ilpm) .860 



= (10.43) (3100) 
= 10,748P‘^^-^ 



Rotor Drive Shaft 



,^,x . u o 7 7 

5.56 h-= = 5. 56(3100) • ^P 

(ilpm V.J.) [j] 

= 1545P*^ 



Tail Rotor 



32 .22 



W 



1.14 



(ilpm V.p) 



1.14 



307,600 

V 

Vt 
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B5.7 Tail Rotor Gear Box 

,.75 

3.7 ^ = (3. 7) (3100r^P' ^ 

(2,pm V^) [^] 

= 206P’^ 



B3.8 Tail Rotor Drive Shaft 



124 



W 



1.355 



(£pm) 



.57 ,,.785 
w 



= (.124) (3100) /A 



A 



. 785 



12.12P-^^ /A 



B3.9 Landing Gear 



Q1 6 99 

191W + .0294W 



B3.10 Rotor Blades Articulated 



19 . 77 



,,1.206 .33 
W a 



DL 



.205 



, , .205 .33 

= 19.77 A “ 

Vt 



B3.11 Rotor Hub Articulated 



.00975 



W 

DL 



1.21 



21 



= .00975WA 



.21 
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B5.12 Fuel System .0615 W 

F 

Calculation of fuel weight three hours at 
cruise SHP 
1515 lbs + 101 
1664 lbs 

B3.13 Total Equation 



WB 



= 12,987,* + 



307600 



V. 



1.14 



+ 206P 



107948P'®^^ + 
^ + 12 . 12P' 



1545P' ^ 

/PT 



+ .191W'^^^ + .0294W^^ 

+ 19.77 + .00975WA‘^^ 



B4 HOVER EQUATION 

Following the formulation in Section of Chapter 3, 
the weight equation based on the design mean lift coeffi- 
cient and power required is: 



W 




411 . 51 



DL 

V. 



3/4 

372 



1 + K, 



Vn 



/ur 



+ Kg /DL 




K 



4 



This number was increased from 8987 to 12987 to bring 
the curves together. This reflects a 4000 lb useful load. 
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where ; 



K, = (1 + 1.807S C,„ 

LRo '' 



C 



LRo 



^2 ^16000/950 



^ _ 0.00025929 ,, , 2^ 

^3 q- (1 ^ 1-^O^^C^Ro ) 



_ 553480.0 
^4 

_ 369 5.7 
^5 

B.5 GR.4PHICAL RESULTS 

Figure B.l is an example of equation (3.13) plotted 
against equation (B.4) for a specific design mean lift 
coefficient . 

Figure B.2 illustrates the family of curves obtained 
when the design mean lift coefficient is varied from 
0.3 to 0.7 . 

In Figure B.3 the solution matrix depicted in Figure 
B.2 is narrowed by the constraints placed on the gross 
weight, rotor diameter and aspect ratio. 
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Helicopter Carpet Plots: CLfl=70 
Utility Class 



CM 

CD 




00961 0006L 00981 

SQl SS0J9 
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Disk Loading 

Figure Bl. Helicopter Carpet Plots: = .70 

Utility Class 



Helicopter Carpet Plots 
Utility Class 



a> 




sqi ;m6i0/v\ ssojq 



o> 

c 

ID 

CO 

O 



cn 



Q 
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Figure B2. Helicopter Carpet Plots 
iJtility Class 



Helicopter Carpet Plots 
Utility Class 



e> 




a> 

c 

CD 

O 

-J 

cn 



sqi m6i0/v\ SSOJ0 
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Figure B3. Helicopter Curpet Plots 
Utility Class 



APPENDIX C. CARPET PLOT METHODOLOGY FLOW CHART AND 
EXAMPLE PROGRAMS: 

This section contains a flow chart to help organize 
a carpet analysis and example IBM computer programs to 
produce the data sets and disspla graphs. 
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«ooly 

Doondary 
condt ton« 



(;♦*♦♦♦♦♦♦♦* GRAPHICiii hELICCPTEF DESIGN FBOGRArt 

C#*#***** CARfrT PLCT ^^U^^EEa 1 ♦♦♦**♦•♦ 

^♦♦♦♦♦♦♦* EY AL HANSEN <c^^***«* 

C*»*^**** this FECGrAK IS DESIGNil TO PBOVIDE A WEIGHT VS. 

C^m****** cisK ICACING GHAFH AT SEVERAL DIFFERENT ROTOR 

TIP VELOCITIES FOR A SPECIFIED DESIGN MEAN 
LIFT COEFfEICENT. 

CBSERVATICi; CLASS 

^4^m******* 4«**4:^* 4 44>«« 4 4c4c«c 4:4*««44««*«* 44:*4i4[««#«4c«««4c4e4e4c4c4:4cft««4>4c««« 44 




VARIABLES: 



Q «« * 4 4« 4«:4 3 

C** ** 

c** 

C< 

c* 
c* 

€♦ 

C* 

C* 

C* 

c* 
c* 
c^ 

c < 



k44444444«4c4«44’i 



NOKENCIfiTUBE 



CLB 

VT 

DL 

WE 

FA 



EESIGN KEAN LIFT COEFfICENT 
TIP VELOCITY 
DISK LCADING 

WEIGHT IS CALCULATED FFCM POWER 
CSEFOL LCAD ILUS EMPTY WEIGHT 
fCWER AVAILABLE IN HCRSEPO'WEE 



k444«*»444*«4«4«44c^444] 



EQUATION 



c4c4(4449c«44 
4 4c 4c 44c 

4 44 
44 
44 
44 
44 
44 

44 
44 
44 
4 4 

4 4*44 



L CLR^IL (15CI,W1 (150) - W2 (15C) /W3 (150) -WU(150) ,W5 (150) 
(150) ,WB 2 ( 156 ) /WB3 (IbS) ,WcMl^ 0 ) ;W£5 (^50) 



C 
C 

c 



SEAL 
1 WB1 



INTEGER I 

DEFINE FILES 
CALL FBTCKS 



1 



( » EILEDEF 
*LATA 



*) 



CALL FBTCKS |*FILEDEr *,*U 
1 * t A T A • * A 

BEAD DATA FROM FILE: CHPTX 



DISK 



DISK 



*CR?T35 » , 



* , *CHPT45 ' 



DATA 



10 

c 

c— 

c 



DC 10 1=1,99 
CRFTl DATA A 

HEAD (3.7C) DL (I) ,W 1 (I) , WEI (I) 
CHFT2 DATA A 
BEAD (4,71) W3 II) ,WE3(I) 
CONTINUE 



2 (I) , WE2 (I) 

W4(I) ,WB4(I) ,W5(I) ,WB5(I) 



-- CALL DISSFLA BCCTINES FOE PLOT 



CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 



TSK6 18 
K EDEUF 
RESET 
HWSCAL 
PAGE (1 
GBACE (O.C 
PHYSCR ' ' 

A BEA2D 
FBAKE 
S WISSL 
BASALF 
MIXAIF 
INT AXS 
YTICKS 
XTICKS 
SHDCHR 
HEIGHT 
XNAKE 
Y NAKE 
HEIGHT 
HEADIN . 



( • SCREEN ») 
2.C, 9.5) 

1 . 2 ) 

1 C-C ,0.5) 



(* L/CSTD*) 
(•STAND*) 



CALL MESSAG 



1 



i-lsf’-' 

(’(D)ISK (L) OADING3 ' , 100) 

('(G) BOSS (W) EIGHT ) (LBS) ) 3 ' , 1 00) 

(.:90) 

(’, iH) EIICOPTEB (C)AEPET (E)LCTS: 
CC, 1.0,1) 

(C)iBFET )F) LOTS: 



015, 1) 



) ' tP) EIICCFTE5 

1CL,3.2E,7. 55) 



(CL3 = 0.5) $' 
(CLB=0. 5) 3’ 
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nonnn non 



70 

71 



CALL HEIGHT (-20) 
k.JxlL ohAH 

1 1 V 3 ALI) = !> 

CALL THKCBV (-C18) 
CALL LEGIIN 
CALL CURVE 
CALL CURVE 
CALL CURVE 
CALL CURVE 
CALL CURVE 
CALL DASH 
CALL CURVE 
CALL CURVE 
CALL CURVE 
CALL CURVE 
CALL CURVE 
CALL THKCBV 



,5^, .2ddu 



(DI,V1 ,95 ,0) 
DL,«2,95,0S 
(Dl,ii3,95,0) 
(DL,'.i4,95,0i 
(DI,Vi5 ,95 ,0) 

(D1 ,1iE1 ,95,0) 
?DL,W£2, 95,0) 
DI ,«E3,95,0) 
DL, W5U,95,0) 



(D1,^E5,95,0) 

(. C30) 

CALL 3LREC r 4 • 2 ,4.6 5 , U 6 , 1. 3 5 , 1) 
CALL GHIC (1, lYGRIC) 

CALL BLOFF (1) 

£3AXLIN = LINEST ( IB AK 1 , 300 , 2 G) 

CALL HEIGHT 
CALL LINES? 

CALL LINES 
CALL LINES 
CALL LINES 
CALL LINES 
CALL HYLEGN 
CALL LEGEND 
CALL DONEEL 
STOP 

FORMAT STATEMENTS 
FORMAT 
FOBMAT 
END 



<2,0l 

* Tfc) EIGHTS* .IRAKI .1) 

* (I) NEUCEDT^ , IPAK 1,2) 

* IF BCFIIES • , IPAK 1,3) 

* (F) ARASITES' , IRAKI, 4) 

( * IP) CWEP (C) URVES* ,1o) 
(IFAK1,4,4.43 ,4.8) 



(5 |2X, F 10.3 
i6|2X,F10 ■ 






76 



^*******4*m***^*3i$**^ ****** 3 

GEAEHICAI hELICCPTEi DESIGN P50GRA.*! 

FAMILY OF SOiailCNS 
CAFFZT PLOT NUMBER 2 
EY AL HANSEN 

THIS P5CG5AM IS DESIGNEI TO ILLUSTRATE TdE Y kli 
CF SCICTICNS FOR HOVEF AND G5EFUL LOAD 
REQUIREMENTS OF A TEETERING ROTOR SYSTE.*- 
WITHCIT ANY BOUNDARY CONDITIONS 
CaSERVATICN CLASS 




■. » * Jtt * ** 

t ****** ** 
********** 
********* 
******** 
******** 
******** 



ILY 



VARIABLES: 



Q** * ** * * **** : 

c** ♦♦ 

C** 

C* 

c* 
c* 
c^ 
c* 
c* 

c* 

c* 



NOMENCIATURE 



CLB 

VT 

DL 

i 

BE 

PA 



EiSIGN MEAN LIFT COEFEICENT 
TIP VELOCITY 
DISK ICADING 

BRIGHT AS CALCULATED EFCM POWER EQUATION 
CSEFUL LOAD ELOS EMPTY HEIGHT 
PCBEE AVAILABLE IN HCESEPOHER 



* * t * ** * * 
******** 
* ** * ^ * ** ** 

t********* 
***** 
*** 
** 
** 
** 
** 
** 

** 
** 
** 
** 

*** ** 



C* 

C* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 



C13 EQUALS THE DISK LOADING FCf CLR= . 3 

C14 EQUALS THE DISK LOADING FOR CLR=.4 

DL5 EQUALS THE DISK LOADING FCR CLR= • 5 

DL6 EQUALS TfcE DISK LOADING FCR CLR=.0 

CL7 EQUALS THE DISK LOADING FCR CLR= . 7 

V3 EQUALS THE WEIGHT FOR CLR=-3 
B4 EQUALS THE WEIGHT FOB CLB=.4 
B5 EQUALS THE WEIGHT FOR CLR=.5 
Wb EQUALS THE WEIGHT FOR CLR=.6 
W7 EQUALS THE WEIGHT FOR CLR=.7 
BVT1 EQUALS nEIGHTS AT VT=600 

■VT2 EQUALS WEJGHTS AT VT=625 

WVT3 EQUALS WEIGHTS AT VT=650 

BVT4 EQUALS WEIGHTS AT VT=675 



c* 


DVT1 


EQUALS 


IRE 


CORRESPONDING 


DISK 


LOADING 


AT 


VT=6U0 


c* 


DVT2 


EQUALS 


THE 


CCHRESPONDI NG 


CISK 


LCADING 


A r 


v;=625 


c* 


DVT3 


EQO ALS 


TfcE 


CCRRESPONDING 


DISK 


LOADING 


AT 


VT=650 


c» 


D VT4 


EQUALS 


THE 


CCRRESPONCING 


CISK 


LOADING 


AT 


7T=675 


c* 


DVI5 


EqU ALS 


TEE 


CCRP.ESPONDING 


DISK 


LCADING 


AT 


Vr=700 



Q*************** **********: 



\ ^ ******** i 



c 

c- 



REAL *4 DL3 (5) .DL4 (5) ,DL5 (5) .Die (5) ,DL7 (5) ,WJ (5) ^ W4 
♦W5 (5) 'W6 (5) (5) ,wVt1(5) ,Ut2 (5) , «V T3 (5) , W VTu (5) , 

^WVT5 (b) ,CVTl (5) ,DVT2 (5) , DVT 3 (5) , DVT4 (5) ,DVT5 (5) 



(5) 



DEFINE DATA 

DATA DL3/2. 06, 2.07,2.07, 2. 07, 2.05/ 

DAIA W 3/^505. , 24 7 6 . 2 , 24 55 . 4 6 , 2433 . 85 ,2418.25/ 

DATA DL4/2.42 , 2.4 4, 2.4b, 2 .4 7 , 2.4 7/ 

DATA S 4/2467. 1 7 , 2U 4 2 . 95 , 2 41 9 . 62 , 2 3 99 . 7 2 , 2 382 . 9 9/ 
DATA DL5/ 2. 6 3. 2. 6 5, 2. 67, 2 .69,2.69/ 

DATA • 6/2440. 72,24 19.51 ,2 399.00,2379 .13,23 64.65/ 
DATA DL6/2. 74 ,2. 77, 2. 79, 2 . 8 1,2. el/ 

DATA Wo/ 24 24. 68, 2402.48, 2 38 3. 5c, 2 36 5. 10, 2 3 52./ 




DATA WVT 4/24 3 3. 6 5, 23 99. 7 2,2 3 79. 13,2365. 1,2357. 35/ 
DATA WVT 5/2 4 1 6. 25 ,23 62. 99,2 36 4.65,23 5 2. ,2 3 42.2 5/ 
DATA DVT1/2.06,2.42, 



3,2.03,2.74 ,2.80/ 

DATA D VT 2/2. 07, i. 44, 2. 65, 2. 77, z. 8 3/ 



C 

C 



DVT 3/2. 0 7 ,2.46,2.67, 2. 79 ,2.85/ 

DVT4/2.07,2.U7,2.69,2.81 ,2.86/ 

DAIA DVT5/2. 0 5 ,2.47,2.69,2. 8 1 ,2. 87/ 



DAIA 

DAIA 
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nnoonnnonnnn nn 



■— CALL DISSELA ECDIINES FOB PLOT 



CALL TEK618 
CALL aZDEUF 
CALL RESfl (3HALL) 

CALL H»SCAL (‘SCREENM 
CALL PAGE (12. C, S.5) 
CALL GFACE (O.Cf 



:0. C) 

( 1 . 0 , 1 . 2 ) 
(1C. 5^0. 5) 



CALL PHYSCa 
CALL AREA2D 
CALL F5AF.E 
CALL SWISSL 
CALL 3ASALF ( » L/CSTE M 
CALL aiXALF ( * SIAAE* ) 
CALL INTAXS 
CALL YTICKS 
CALL XTICKS 
CALL SHDCHB 
CALL HEIGET ( . 1fa[ 



hio, 1,. 



015,1) 



CALL XNAHE (*(E)isK (L) OA DI N G 3 » , 1 00) 
CALL y SAKE »(g[eCSS (W) EIGHT ((LBS 
CALL HEIGHT (.290) 

^CALL HEAEIN ( »^ ^H) ELICCPT E fi (C) ARPET 
^CALL HEIGET (.^20)" 

CALL GRAF ( 2. 0 , . 1 , 3. C , 2 3 0 C. , 5 C . , 2 550 
iyGRID=5 

CALL THKCBY (.CIS) 

CALL PABA3 
CALL LEGLIN 
CALL CURVE 
CALL CUBVF i 

CALL CURVE i 

CALL CUaVE I 

CALL CUBVE 
CALL DASE 

CALL CURVE (D V1 1 , 6 VT 1 ,5 , 0 ) 

CUBVE (DVT2,iVT2,5, 0) 



DL3,Si3,5 ,0) 

DL5,«5,5,0) 
Oi6,ii6,5,0i 
DL7, *7,5,0) 



CALL 

CALL CaaVE (DVT3, iVT3,5 ,0) 

CALL CUaVE (DV14,IVT«,5,0) 

CALL CUaVE (DVl5.iVT5,5,0) 

CALL THKCaV ( . C30) 

CALL BLaEC 14.2,4.65,1.6,1. 35,1) 
CALL GBIC (1, lYGHlC) 

CALL 3LOEF (1) 

MAXLIS = LINEST ( IP A M , 300 , 2 0) 

CALL HEIGHT (. 12) 

CALL LINESP (2.0) 

CALL LINES ( • ( i) EIGHTi' » I PAK 1 , 1) 
CALL LINES ( ’ ( I) N COCEDi • , IP AK 1 ,2 



CALL LINES ( ' ( P BCFI LE5 ' , IP AK 1 , 3) 
CALL LINES (• ) E 1 A r AS ITES ' , IP A K 1 , 4 
CALL SYLEGN r(P)C'«EF (C)l)EVES',1 
CALL LEGEND (I P AK 1 , 4 ,4. 4 3 , " 



CALL DONEEL 
STOP 

C ECdKAT SIATEHENIS 

C70 PCRHAT (5 (2X,F10.3)) 

C71 FOaaAT (6 (2X, F 10.3) j 

END 



. 8 ) 



1^) 



) S' , 100) 

(P) LOTSS • , 

-) 
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«V«:ti««i;r«)ic4c««4e4c4*««4cKtt*444*A:kA4; 

GnAtHICfll hSLICCPTiif DESIGN F 50 GHAM ♦ » 4:4c ^ 44 

ASPECT RATIO BOUNDARY 4 * 44^444 

LOCI Cf HCVER AND USEFUL-LCAD SOLUTIONS 4444*444 

CABF 5 T PLOT NCtcER 3 44444444 

BY AL HANSEN 4*444444 

44444444 

TBIS PBCGEAK IS DESIGNEE TO GRAPHICALLY DETERMINE 44444*44 
THE ASPECT PATIO bOUNEARY REQUIREMENTS FOR 4 * 4444*4 

A HCTCR SYSIEW USING EBEVIOUSLY GENERATED DATA 4 * 44 * 44 * 

444444444 * 

444444444444444444444 = 

44 444 4444 4*4 3 



NOaSNClATURE 



COEFEICENT 



EESIGN BEAN LIFT 
TIP VELOCITY 
DISK ICADING 

ASPECT PATIC. HISTOEICALLLY 



44*44 
4 44 
44 
44 

44 

44 

44 

ASSUMED TO BE LESS THAN 21 

44 
4 4 4 44 



DVT1 

DVT2 

DVT3 

EVT4 

DV15 

C600 

C625 

C650 

C675 

C700 



sguAis 

EQUALS 

EQUALS 

EQUALS 

EQUALS 

EQUALS 

EQUALS 

EQUALS 

EQUALS 

EQUALS 



TEE 
Tl E 
THE 
TEE 
TEE 
TEE 
TEE 
TEE 
TEE 
TEE 



CCSRESPCNDI NG 
CCFRESPONCING 
CCRRESPONCING 
CCREESFONDING 
CCBRESPCNCI NG 
LIFT COEF AT 



LIF1 

LIFT 

LIFT 

ur 



COEF 

COEF 

COEF 

COEF 



AT 

AT 

AT 

A 



EISK 

DISK 

DISK 

DISK 

DISK 

VT=bOO 

VT=d25 

^1=050 

VT=o75 

\T=700 



LOADING 

LOADING 

LCACING 

LOADING 

LOADING 



M X 

AT 

AT 



VT=600 

VT=625 

VT=o50 

VT=675 

VT=700 



144*4 4*44444444 444444 4444 ’ 



(444444444444 



BEAL*a 

♦C675 
♦DV11 



CLF (5) ,EL (10) ,C6 0 0 (10) ,C625 (1 C) ,C650 ( 10) 



n 0) ,C 700 ( 10) f 

(5) /CVT2 (5) ,dCt3 



(5) , D VT4 (5) , DVT5 (5) 



C 
C 

C 



DEFINE DATA 

DATA DVT 1/2. 06 ,2.42,2.63, 2. 74 ,2.80/ 

D VT 2/2. C7 ,2.44,2.65 ,2. 77,2.33/ 

DVTj/ 2. 07,2.4t,^.67, 2. 79,2. 85/ 

DVT 4/2. 07, 2. 47, 2. 69, 2. 81 ,2.3o/ 

DVT5/2. 05 ,2. 4 8, 2.7 0, 2. 32 ,2-8 7/ 

CLR/.3,. 4 ,-5,.o ,.7/ 

DL/2-,2. 1,x.^, 2. 3, 2. 4, 2. 5, 2. 6, 2. 7, 2. 6, 2. 9/ 
G60C/.47 ,-49, .51, .53 »«55 

C62 5/. 43, .45, .47, -49, .51, -53, .55, .57, .59, .61/ 

C 65 C/. 3 9,. M,. 4 3, . 4 5 , . 47 ,. u9 , . 5 1 ,. 5 3 , . 5 5 , . 5 7/ 
C675/. 35,.j7,.j9,.4 1,.43 ,.45,.47,.49,.51,.d3/ 
C70C/. 31, -33, .35, . 37, .39, .41, .43, . 45 , .47, .49/ 



DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 



CALL DISSPLA fiCUTINES FOR PLOT 



CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 



.0,6.5) 



TEK6 18 
MEDEUF 

RESET (3HALL) 
aWSCAL ( * SCREEN •) 
PAGE (1^-0, 9.5) 
GRACE (C.0‘ 

PHYSCR (1- 
ABEA2D (10 
F SAME 
SWISSL 
BASALF 
MIXALF 
INTAXS 
S HDCER 
HEIGHT 
X NAM_ 
y NAME 
HE 

MESSAG 



iL. J.VJ ni { • ' u 

:name (* (D) 

'NAME (• \C\ 
iEIGET (.29 



(•L/CSTDM 
( * STAND* ) 

(.90, 1,- 015, 1) 

(. 16) 

£) ISK (L) OADING3 *, 100) 

' CZEFICENT OF (1) IFTi 

0 ) 

(H) ELICOPTER (C) ABPET 



• . 100 ) 

(?) ICTSS > 
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1 ICO. 2-25, 7. 

CACi. li Lxu Cl (. ^0 I 

L k L±, k At ^ , km k / • k 

iyGRID=5 

CALL THKCBV (.CIS) 

CALL PABA3 
CALL LLGLIN 

CALL CURVE (D I ,C6C C , 10, 0 
CALL CURVE ( D1 , C6 25 , 10, 0 
CALL CURVE ( D1 ,C6 50 , 10, 0 
CALL CURVE { D1 , Co 75 , 10, 0 
CALL CURVE (Di ,C7CC , 10, 0 
CALL DASK 

CALL CURVE (D V T 1 , CL3 ,5 , 0 
CALL CURVE (D V 1 2 , CLS , 5, 0 
CALL CURVE (D VIS , CLR ,5 , 0 
CALL CURVE ( D VTm , CLR ,5 , 0 
CALL CURVE (D V 15 , CLB ,5 , 0 
CALL THKCSV (.C30) 

CALL DCNEEL 

STOP 

EHC 



5 5) 

,-05, .t) 

1 
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C ** <*»♦♦♦ 
Q*:$i t ** * 



c ♦*<<♦♦♦* 
c**<*^**< 

c## ««««## 



► »»♦»» 4 4* 4 4 *♦ 4jp**<t» 4X«)pj«t«9«4c«4(Jt(4i4i4c4(i»«94C9«^« 

GRAPHICAL hiLiCLFTtf DESIGN F50GKA!! 

PAtILY Or SOLUTIONS **44*4r4c4* 

CAEPET PLOT NUKEER 4 
EY AL HANSEN 

THIS PECGEArt IS DESIGNEE TO ILLUSTRATE THE FAHILY 
CF SCIOTICNS FOR HOVEF AND LSEfUL LOAD 
EEQUIEEEEMS OF A TEETERING FOTOR SYSTSi- 
-ITH ECTCE EIANETER AND KAX GROSS 30UNDRIES. 

CBSEBVATION CLASS HELICCFTFfi 

( 4 ♦ 3 



C4#44^44*^#*#4*##*#444*J 
Q t* n t* 1 1 ****** t *** * t * n ^ n n ********** n nti **** ^ 
C** ♦♦ 

NOMENCLATURE 



C** 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c*** ^ 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c^ 

c* 

c* 

c» 

c* 

c* 

c* 

c* 



VARIABLES: 

CLE EESIGN MEAN LIFT COEFEICENT 

VT TIP VELOCITY 

DL DISK LOADING 

m HEIGHT AS CALCULATED EECM POiiER EQUATION 

■B USEFUL LOAD ELDS EMPT^ WEIGHT 

PA POWER AVAILABLE IN HCRSEPOWSB 




ELS 

CI4 

DL5 

CI6 

EL7 

K3 

W4 

W5 

W6 

H7 

W VT1 
WVT2 
W VT3 
tfVTU 
VT5 
EVT1 
DVI2 
DVT 3 
DVT4 
DVT5 
RES 
WKG 



EQUALS 
EQUALS 
EQUALS 
EQUALS 
EQUALS 
EQUALS 
EQUALS 
EQU ALS 
EQU ALS 
E^U ALS 
' EQUALS 
EQUALS 
EQUALS 
EQUALS 
EQUALS 
EQUALS 
EQUALS 
EQU ALS 
EQUALS 
EQUALS 



THE 

ThE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

HE 



DISK 
DISK 
DISK 
DISK 
DISK 
WEIGHT 
WEIGHT 
WEIGHT 
WEIGHT 
WEIGHT 
WEIGHTS 
WEIGHTS 
WEIGHTS 
WEIGHTS 
WEIGHTS 
Tfc 
EE 



FOE 

FOR 

FOE 

F.CE 

FOE 



LOAD I NG 
LOADING 
lOADI NG 
L 0 A D I NG 
LOADING 
FOR CLR = .3 
FOR CLR=.4 
FOR CLR=.f 
FOR CLR=.6 
FOR 01? = .“/ 
AT VT=6 0 0 
VT=o25 
VT=650 
VT=o75 
VT=700 



CLR=.3 
CLR=.4 
CLR=. 5 
CLR=.o 
CLR= .7 



AT 

AT 

AT 

AT 



TEE 

TEE 

'EE 



DISK 

DISK 

DISK 

DISK 

DISK 



LOADING 

LOADING 

LOADING 

LOADING 

LOADING 



C 

c 



CORRESPONDING 
CORRESPONDING 
CORRESPONDING 
CCRRESPONCI NG 
CCHR ESPONDING 
ROIOR DISK ECONDAHY 
MAX GROSS WEIGHT BOUNDARY 

1t*:******H 4 *if*** *** 

BEAL^4 
♦W5(5) 

♦ W VID 



AT 

AT 

AT 

AT 

AT 



VT=600 

VT=d25 

VT=d50 

VT=d75 

VT=700 



♦4 DL3 (5) ^DL4 |5| ,DL5 (5) - DL6 (5) ,D L7 (5) , W3 (5) . W4 (5) , ROB 
. Wo (5) , « / (5) ,W VT1 (5) ,WVT^ (D) ,WVT3 (5) , WVtu (5) ,WMG (2) , 
(b) ,CVT 1 (5) ,DVT 2 (5) , DVr3 ( 5) , DVT4 ( 5) , D VT5 (5 J , D 1 ( 2) , D2 ( 2 ) 



33(2) 



DEFINE DATA 
DATA DL3/2.06 



DATA DLJ/^i.Ub^i'.L / , 2. 07 , 2. 07, 2. 05/ 

DATA W 3/2505. , 2478. 2 ,24 55.46,2433.85 ,24 18.25/ 
DATA DL4/2. 42, 2. 44, 2. 40, 2. 47, 2-47/ 

DATA W 4/ 2 467. 1 7,244 2.95, 24 19. 62,2399 .72 ,2382.9 9/ 
r a'T’ a n t ^ ^ Z ^ 1 1 Q q / 



1ST A DLd/ 2. 63,2.65,2.67,2.09,2.09/ 

DATA W 5/ 2 440. 7 2, 24 19. 51 ,2399.00,2379 . 13,2364.65/ 
DATA DL6/ 2. 74-2. 77, 2. 79, 2. 81, 2. 01/ 

DATA W 6/2 4 24. 0 6, 2402. 48, 2383. 56,2365 . 10,2352./ 
DATA DL7/2.80,2.8^,2.65,2.8o,2.o7/ 

DATA W 7/ 241 1. 92,2390.72 , 2372. 50,2357 .35,2342.25/ 
DATA WVT 1/250 5. ,246 7. 17, 2 44 0 . 7 2 , 24 24 . 67,24 11.92/ 
DATA WVT 2/24 7 8.2, 24 4 2.9 5 ,24 1 9 . 5 1 ,24 0 2 . 5 3 , 2 39 0 . 72 / 
DATA WVT3/24 5 5. 4o ,24 19. 6^,23 9 9., 2383 . 56 ,2372. 5/ 
DATA WVT4/24 3 3. 6 5 ,2399.72 ,2379. 13,2365. 1,2357. 35/ 
CiTA HVT5/2ine.:5,23e2. 99, 2 36U. 65,23 5 2. 342.25/ 

DATA DVT 1/2- 0 6, 2. 42,2- 63, 2. 74 ,2.80/ 

DATA DVT2/2- 07, 2. 44 , 2-d5 , 2- 77 ,2.83/ 

DATA D VT3/2. 07 ,2.46, 2.6 7 , 2. 79 ,2.85/ 

DATA DVT4/2. 07 ,2. 47, 2.69 , 2. 61 ,2. 86/ 

DATA DVT 5/2. 0 5,2.47,2.69,2.81 ,2.87/ 
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onoo 



DATA 01/2.52-2.?/^ 
Lili-A wrili/^4D0* f 2 ^ 5 u • / 
DATA D2/2-J6 - 2.-2/ 
DATA RDB/230U. , 2*4t0. / 



CALL DISSELA BCDTIl.ES FOR PLOT 



;iC. ,6.5) 



I JI/CSTC *) 



»£TA^D * ; 



(.50,1,. 



015,1) 



CALL TSK618 
CALL KEDEUF 
CALL RESET (3HALL1 
CALL HWSCAL MSCRrEN*) 

CALL PAGE (12. C, S.5) 

CALL GRACE (0-0)i 
CALL PHYSCE (1.0,1.: 

CALL A REA2D 
CALL FRAKE 
CALL SWISSL 
CALL BASALF 
CALL MIXAIF 
CALL INTAXS 
CALL YTICKS 
CALL XTICKS 
CALL SKDCKR 
CALL HEIGHT j. 16) 

CALL XNAHE (*(L)ISK (L) 0 A DI N G E • , 1 00) 
CALL YNAHt (*(G)RCSS (W) EIGHT ( (LBS) 
CALL HEIGHT (. 250) 

CALL HEADIN (M H ) ELI COPT S H (C)ABPET 
1 ICC, 1.0,1) 

CALL HEIGHT (.20) 

CALL GRAF (2. C , . 1 ,3. C,2 3 0 C. , 5 C . , 2550 
IYGBID=5 

CALL THKCEV (.C 18) 

CALL PARA3 
CALL LEGIIN 

CALL CURVE (DL3, h3,5 ,0) 

CALL CURVE (DL4,«m,5,0) 

CALL CURVE (DL5, '..5,5 ,0) 

CALL CURVE {DI6,Wc,5,0) 

CALL CURVE ( D I V , Vi V , 5 , 0) 

CALL DASH 

CALL CURVE { D V T 1 , V. V T 1 , 5 , 0 ) 

CALL CURVE { D V T 2 , h V T 2 , 5 , 0 ) 

CALL CURVE ( D V T 3 , V. V T 3 ,5 , 0 ) 

CALL CURVE ( D V T*^ , V V T 4 , 5 , 0) 

CALL CURVE (D VT5 , ; VT5,5 , 0) 

CALL THKCEV (. C30) 

CALL CHNCCT 

CALL CURVE (D1,WKG,2,0) 

CALL CURVE (D2, RLE, 2,0) 

CALL DONEEL 

STOP 

ENC 



) E* ,100) 

(r) LOTS5 * , 

.) 
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q4#4 4e#««4 

Q^m ♦♦♦♦*♦ 



c* 

C» 

C* 

C* 

c* 

c* 

c^ 

c* 

c* 

c* ♦♦♦ 

c» 

c» 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 



VARIABLES: 



. mm mm mm #»»»«♦# ^ 4 - 4 

GRAPHICAL hELICCPTLr bcjlo « ttUbwAfl ^ m- ^ m - 

rAEILY OF SOLUIICNS »,444444« 

CABEET ?LOI NU^lcER 5 ««44**** 

£Y AL HANS£^ :i« « ♦ :* <t 

4 C 44444 C*# 

THIS PRCGBAM IS DESIGNcC TC ILLUSTRATE THE FAMILY «♦««##♦« 
Cr SCIUTICKS FOR HOVER AND LSEFUL LOAD ♦♦♦♦♦«♦♦ 

REQUIREMENTS OF A TEETERING ROTOR SYSTEM 
WITH RCTCR DIAMETER, MAX GROSS WEIGHT AND 
ASPECT RATIO SCUNDARIES. 

CESERVATION CLAES HELICOPTER 

4 ♦ ♦ ♦ « 

NOMENCLATURE 

♦ ♦ 

4 ^ 
44 

4* 




CIB DESIGN MEAN LIFT COEFilCENT 

VI TIP VEICCITY 

Cl DISK LCACING 

W WEIGHT AS CALCULATED FfCM POWER EQUATION 

WE OSSFUL LOAD PLUS EMPT^ WEIGHT 

PA PCWER AVAILAcLE IN HORSEPOWER 



CL3 EQUALS THE DISK LOADI-i^G PC? CLH=.3 

DL4 EQUALS THE DISK LOADING FOP CLR=.4 

EI5 EQUALS THE DISK LOADING FOR CLR=.5 

DL6 EQUALS THE DISK LOADING FOB CLR=.6 

EL7 EQUALS THE DISK LOADING FOB CLR=.7 

W3 EQUALS THE WEIGHT FOE CLR=.3 
WU EQUALS THE WEIGHT FOR CLB=.4 
W5 EQUALS THE WEIGHT FOR CLR=.5 
W6 EQUALS THE WEIGHT FOP. CLR=.6 
W7 EQUALS THE WEIGHT FOR CLR=.7 
HVTl EQUALS WEIGHTS AT VT=600 

WVT2 EQUALS WEIGHTS AT VT=625 

WVT3 EQUALS WEIGHTS AT VT=650 

WVT4 EQUALS WEIGHTS AT VT=675 

WVT5 EQUALS WEIGHTS AT VT=700 



44 C 

44 
♦ « 

4 4: 

4 4 « 44; 



C* 


DVTl 


EQUALS 


THE 


CCRRESPONDI NG 


DISK 


LCADING 


AT 


VT=600 


c* 


DVT2 


EQU ALS 


Tt E 


CCRR ESPONDING 


DISK 


LCADING 


A T 


VT=b25 


c* 


DVT3 


EQUALS 


THE 


CCERESPONDING 


DISK 


LCACING 


AT 


YT = d5 0 


c* 


LVT4 


EQU ALS 


THE 


CCR5 ESPONDING 


DISK 


L C A D I N G 


AT 


VT=o75 


c* 


DVT5 


EQUALS 


THE 


CCRRESPONDI NG 


DISK 


LCADING 


AT 


VT=700 


C* 


RDE 


ROTOB DISK 


ECUNCARY 











c* 

c* 

c* 

C444 ; 



c 

c 



WMG MAX GROSS WEIGHT BOUNDARY 
AR WEIGHTS FOR THE ASPECT RATIC BOUNDARY 

D3 DISK LOADING CCC B ESP 0 N DI NG TC IKE ASPECT RATIC WEIGHTS 

44444»4»4t44444 4 44 4 444 4*4:4 444 4*4 44*444 444444444*44 4444444»44 44 444444 

RE&L*4 DL3(5).CLU(5),DI.5(5),Dltf5),DL7(5),V.3(5).»U(5) ,HDE (2) , 
*»5J5) • '<6 (5i,ii7 (5) .iVII (5/ ,«'VT2 iS) ,WVT3 (5) ,SVT4 (D) .aKG (2) , 

*HVI5 (5) ,EVT1 (5) ,fcvi2 (5) , D VI 3 ( 5) , DV74 (5) ,d0t5 (5) ,D1 (2) ,D2 (2) , 

*D3 (5) , AR (5) 

DEFINE DAIA 

DATA DL3/ 2. 06, 2. 07, 2. 07, 2. 07, 2. 05/ 

DAIA U3/2505. , 247£. 2 ,24 d5 .4 6 , 2 433. 85 , 24 18 . 25/ 

DATA DL4/ 2. 42, 2. 44, 2. 46, 2. 47, 2. 4 7/ 

DAIA W 4/2 467. 17,244 2.95 ,241 9.62,2399 .72,2382.99/ 

DAIA 0L5/2.63,2.65,2.67,2.o9, 2.0 9/ 

DAIA W5/2 440. 72,2419.51 , 239 9 . CC ,2379 . 1 3 , 2 3o4 . 65/ 

DAIA DL 6/ 2. 74, 2. 77, 2. 79, 2. 8 1,2. 81/ 

DAIA W 6/2424.06, 24 C2. 4 3, 2333. 58,2365 .10 ,2352./ 

DAIA DL7/2-80 , 2. 6 3 , 2 . 65 , 2. 8o , 2. 87/ 

DAIA S7/24 11. S2,2390 .72 , 237 2. 5C , 2357 . 35 , 2 342 . 25/ 

DATA iiVT 1/2505. ,246 7. 17, 2 440. 72 ,2424 .67,24 11 . 92/ 

DAIA W VI 2/24 78. 2, 244 2.9 5 ,24 1 9 . 5 1 ,240 2. 5 3, 239 0. 72/ 

EAIA W7I3/2455.40, 24 19.0/, 2399. ,2383. 58, 2372-5/ 

DAIA V. VI 4/ 24 3 3. 6 5, 23 99. 7 2 ,2379.13,2365. 1,2357. 35/ 

DATA WVT5/241 6- 25, 23 62. 99,2 364.65^23 52. ,2342.25/ 

DAIA D VI 1/2 -06, 2-42, 2. 63, 2. 74,2.8 0/ 
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nono 



83/ 

CD/ 



DAIA DVT2/2. 07, 2.^44,2.65,2. 77,2- 
l2Ta L>v'ij/z.Li7,ii-4t,2.67, z. 
bit! A Dva4/^.U / ,/.47,2.6':y,Z.t3 I ,Z. dO/ 

CAXA DVT5/2. 05,2.47,^.69, 2. 81 ,2.87/ 

DAIA DI/2.52.2.9/ 

DATA y MG/2450. ,2450./ 

DAIA D2/2.36 , 2.52/ 

DATA aDB/2300. ,2450./ 

DATA A H/ 2 374. 461, 23 9 G., 24 C5., 24 20., 2 450./ 
DAIA D3/ 2. 44, 2.65,2.71?, 2.745 ,2.78/ 



CALL DISSDLA FCCTINES FOR PLOT 

CALL TEK618 
CALL MLCEDF 
CALL RESET (3H£LL) 

CALL HWSCAL ( 'SCREEN •) 

CALL PAGE (12. C, 9.5) 



CALL GHACE ’ (0. C[ 

CALL PHYSCR ( 1.0 , 1. 2) 
CALL ABSA2D (1C. ,6.5) 



{ii 



CALL FRAME 
CALL S»I£SL 
CALL BASALF ('L/CSTC*' 

CALL MIXALF ( • £TA^C *) 

CALL INTAXS 
CALL XTICKS 
CALL XTICKS 

CALL SHDCHE ( . 90 , 1 , . C 15 , 1 ) 

CALL HEIGHT (. 1b) 

CALL XNAME (»(D)I£K (L) 0 A DI N G 3 • , 1 00) 

CALL YNAKE ('(GliBCSS (.4) EIGHT | (IBS) ) 
CALL HEIGHT (.290) 

CALL HEALIN’ ( » JH ) ELI COPT E 5 (C)AB 
1 1CC,1.C,1) 

CALL HEIGHT (.20) 

CALL GRAF (2 . C , . 1 , 3. C, 2 3 0 C. , 5 C . , 2550 . ) 
IYGRID-5 

CALL THKC3V (.C13) 

CALL PARA3 
CALL LEGLIN 

CALL CURVE (Dl3,y 2,5 ,0) 

CALL CURVE ( D L 4 , 4 , 5 , 0) 

CALL CURVE (D L 5 , 5 , 5 , 0) 

CALL CURVE (DL6,Wc,5 ,0 
CALL CURVE jDi7,W7,5,0) 

CALL DASH 

CALL CURVE (D \ T 1 , is VT 1 , 5 , 0 ) 

CALL CURVE ( D V T 2 , V V T 2 , 5 , 0 
CALL CURVE ( D \ T3 , > VI 3 ,5 , 0 ) 

CALL CURVE { D VT 4 , y VT 4 , 5 , 0) 

CALL CURVE (D V T5 , VT5,5 , 0) 

CALL THKCRV (.030) 

CALL CHNDCT 

CALL CURVE (D 1 , WKG, 2,0) 

CALL CURVE (D2, REE, 2.0) 

CALL CURVE (D3,AE,?,C) 

CALL DONEPL 

STCP 

END 



, 100 ) 

(?) LOTSB* , 
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no nno non no no nnnnn n non no nn onnnnr. 



C THIS PROGHAK IS CISIGN TO GENERATE THE lATA FOR THE GRAPHTCAT 
L a»-i.UiiON Of i-rc. AN'» THF USEFUL LCAD EOUATlUN. tmt'- tc 

rinil SXEt' A oAbirFT PLOT HELICCFT2H CESIGN PARAMETRIC QPTIMI2ATICN 

ASSUMPTIONS: 

1> FNGINFSPECIFIEr 

VABIAELE OPTIONS 

BEAL*4 C1E,PA ,EI K3 ,K4 , K5 ,R,S, A,P ( 10) , WB ( 10) 

INTEGER VT,D,i,Cl 



CALL FRTCMS ( • EliECEE ' 
>*CATA *,*A •) 


' , * 02 


• DISK 


*CHPT 1 * 


CALL FBTCKS ( 'FILECEE ' 


1 1 
1 1 
1 1 

1 IS 

1 1 - 

1 1 


* , • DISK 


* , 'CEFTP • 



>*EATA *,»A *) 



CLfi= DESIGN MEAN LIFT COEFFICIENT 
DC 90 CL=3,7 
CLfi=CL^ (C- 1) 

8FITE (2, 10) CLE 

FA= POWER AVAIIAEIE EP 

PA=205 

£1= DISK LOADING 

VT= TIP VELOCITY ET/SEC 

CONSTANTS BASED CN CIE 

K1= (0. 95e31/CIE» ( U 1 .607 8 ♦CLS^*2) 

K2=PA*10**b/K 1 

K3= f 0. 00C2592 9) /CIR» (U 1 . 807 8 CLR**2 ) 

KU=55348C,0/K 1 
K5=3695. 7/K1 
DC 100 E=200,300 

LL-D* (0.C1) 

1 = 0 

DO lie VT = 6CC ,70C,25 

ARRAY INCREMENTER 
1=1 + 1 

WEIGHT EQUATION 

W (I) = (K2* ( 1- (4 1 1. 51 ♦£!** . 75) / (\T^*1 . 5) ^ ( 1 + K3 * VT/DL* 5 ) »♦. 5) -N4) 
1/ (VT + K5*C1^*. 5) 

CALCULATICN OE WE DATA 



P=A*^. 5/VT 

S= (6. ♦ CL)/ (0. 0 02367 9«CLR *VT**2) 
ASSUMING A TEETERING SYSTEM 



WE (I) =1717. 9 + 1 : 21 . ♦P^ + .ae 3+ 266. *?*♦. 7+10449. / VT* *1.14 
1+2,886*A**.5*P**.57+.191*W(I) +*.916+.0294*«(I) **.99 

.... ... ..... .......... .. (I)*A’« 



2+35-15* 



(W(I)/7T)*U**.185)*S**.33 + .0C83*W(J 
WRITE ib, 2C) VT 



.21 



110 



100 

90 

31 



HRITI (5 ,i0) S 
CONTINUE 

(2,3 1 
3,31 






WRITE 
WRITE 
CONTINUE 
CONTINUE 
FORMAT STATEMENTS 
EOEMAT ( ' 1« , • CLfi= • , IE 10. 4///) 
FORMAT (6 (21, E 10.3) ) 

STOP 

END 



E (2) 

"E(4) ,W(5) ,WB(5) 
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C XL2*^ cAas^ounu xo LilDavrti i.o oi:IAx» inZ oAuS::) HiilGhr f hUO TnL ^ 

C ASPECT HATIC GSAFK INTEPSECTION tCINT DATA. IN ORDER TO GRAPH ♦ 

C THE ASPECT RATIO ECUKCAFY ON THE KAIN CARPET PLOT A ♦ 

C CCHEESFCNCING WETGET MIST BE OBTAINED. THIS FROGRAL UTILIZES 

C THE ^iEIGHT ESTIMATICN ECHMOLATIOh FROM THE CARPET PROGRAM. ♦ 

Qif^ if *:$**^if t 4 t * ^ t t t t t n^*^Hf**^***--in**i^* rn^****** 

C ASSUfiPTIONS: ♦ 

C VAfilAELE OPTIONS 

REAL+4 VT <5) .C1J5) ,CIR(5) ,W(5) 

CAT A CLR/.625 ,.5‘}5,.5 35, .4b0,.y2/ 

DATA VT/6CD. ,625.,6 5C. ,67 5. ,7CC./ 

DATA DL/ 2. 75, 2.75 ,2. 72, 2. 67, 2.53/ 

CALL FBTCKS ( » EILEDEF 
>*DATA »,»A ») 



' 03 



DISK 



•CHPTAR* , 



C 

C 

C 

c 

c 



DC 90 L= 1.5 
SRITE(3, 1C)CLP (I) 

^J£ITE(3, 2C VT (I) 

^’BITE (3 ,30) DL )l) 

PA= POVEE AVAILABLE HP 
PA=206 

CL= DISK LOADING 

VT= TIP VELOCITY FT/SEC 

CONSTANTS BASED Ch CIB 

K 1= (0. 9563) /CL5 a) ♦ ( 1+1 . 807 8*CLH (L) ♦ *2) 
K2=PA» 10^^b/K 1 

K3= (0. 00C259 2 9) /CIR (L) « ( 1 +1. 8 C78*CLR (L) *<‘2) 
K4-553UBC.0/K 1 
K5=3695. 7/K1 

HEIGHT £2UATIC^ 






= (K2* (1- 



(yi1.51*CL(i 
CL (1) • * ' 





'-BITS 


(3,31) W (I) 


90 


CCNTINI 


JE 




FORMAT 


ST ATEMENTS 


10 


FORMAT 


(• 1 • •CIR = 
2X, ^ VT = * , 


20 


FORMAT 


30 


FORMAT 


(2i,*DI=* , 


31 


FORMAT 


(2i,F1 C.3) 



•».75)/ (VT ( 1) 

) *■* 1 5) - KU)/ (VT ( L) + K5* (DL (L) **. 5) ) 



= *1.5) 



STOP 

ENC 
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APPENDIX D. PROGR.AMS TO ACCESS HESCOMP 



This section 
needed to access 



contains the control language 
HESCOMP on the IBM main- frame 



programs 
computer . 
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•■ ♦ 


♦ 






OzT 










®3 














• 






• *-uo 




• 


— *^mzr-3 . 


♦ 


♦ 


* 


Oz ♦ ♦ 


♦ 




• o 




*- r-iTunr* 


O !7'C‘ ®004^-n® 




. o 


®mpk^ 3— >30 


♦ 


* 


♦ 


6-^ ♦♦ 


4fk 




o o 


^XZ*-^^<NO fNfNfN*-*-Z 


o xo4trr*c^oz*-^ 




o o 


-OXX— 30k0JO 


♦ 


» 




Z ♦ 




• 


.X . 


.fN^rr- • . • .o . « 


lO— oj-ncrmor- 


® 


• 


m .CO . . 


X3r- * . • .® 


* 


♦ 




Wz ♦ 


♦ 


o 


:yvZ»-0 z 




*- m o o o o o *-•-*“ «- 


o 


!y^O — o o 




♦ 


♦ 




>•3 ♦ 


* 


















♦ 


♦ 




♦ 


♦ 


















♦ 


* 




= * 


♦ 


















* 


♦ 




Cm^ * 


* 


















* 


« 




* 


♦ 


















♦ 


♦ 




O ♦ 


♦ 


















♦ 


♦ 




zto ♦ 


*■ 


















♦ 


♦ 




&2 ♦ 


♦ 


















♦ 


« 




Oz ♦ 


« 






irx 












♦ 


« 




Om ♦ 


« 


=r 


fN « • 


^-NO » • 


• 


« • X iT 




04 • • 


— 3— 04040 . 


♦ 


•* 




-3 « 


♦ 




*^oo 


;»jy'^srmo3’oo 


r-oD®m^o T>o 


OO XX'^®®tOO'^ 




'*noo 


30 — 3 mo 30 


♦ 


♦ 




- •» 


♦ 


\ 


ozz z 


O*-ZX'NJ^^OO'.0OfNOC^ XZ-0*“OZ 


oo® ojino*zoo*“»- 




OOO z 


O— Z o 'Nm3 o 


♦ 


* 




i-TE- * 


« 


^ • 


«oo • « « 


i-'N^OC' « . t^TO « 


iO’-oi'N'ocrinr*' o ♦ 


• • 


• 


•ox . * » 


04303 * • *3 


♦ 


* 




uz « 


♦ 


CN(N 


OX*-UlOO 


• t • •*-*-»-*-rsio 
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APPENDIX E: HESCOMP INPUT AND OUTPUT EXAMPLES 



This section contains samples of the IBM computer 
input and output. 
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